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1. Introduction 

Measurements of the flavor-changing, neutral current (FCNC) decay b —y s are rapidly growing 
in number. In particular, the semi-leptonic FCNC decays of B and B s mesons to the vector mesons, 
K* and (ft respectively, have afforded us the opportunity to compare several associated observables 
with theoretical predictions [3, 4, 5, 6]. The experimental data presently show a few deviations from 
the Standard Model, including a smaller-than-predicted differential branching fraction in low-recoil 
kinematic bins. These could be explained by beyond-the-Standard-Model (BSM) contributions to 
one or two Wilson coefficients in the effective b —> s Hamiltonian [7, 8, 9, 10, 2, 11, 12]. However, 
none of the analyses claim the Standard Model is an unacceptable description of the data. In fact it 
is not yet certain that all theoretical uncertainties are under firm control. 

The largest theoretical uncertainties in these observables are due to QCD interactions. Lattice 
QCD (LQCD) calculations of the Br s \ —y K*/(ft form factors can reduce a class of these uncertain¬ 
ties, improving upon the determinations from sum rules. We recently completed a necessary first 
step toward accurate, first-principles calculations of the form factors [1]. Below we briefly summa¬ 
rize the results of computing matrix elements of the full basis of b —> s' currents. The key improve¬ 
ments made compared to previous lattice studies come through using high statistics, physical-mass 
bottom quarks, and 2+1 flavors of sea quarks. In Sec. 3 we provide a minor improvement to the 
method used to fit the form factor shapes. We close in Sec. 4 by touching upon issues which must 
be confronted in order to improve the theoretical determinations of these observables. 


2. Summary of the matrix element calculation 


The effective Hamiltonian governing b —y s decay is 


'T/sib—ts 

^eff 


^V*V tb 'E(C i Oi + C i O'J 


( 2 . 1 ) 


For the following operators in matrix elements factorize into local hadronic matrix ele¬ 

ments, which are parametrized by seven form factors: 




(?(') = 




sy^Pu^bly,, 


where P L / R = |(1 + y 5 ) and cF 11 ' = ft[y- u .y v . The bulk of this write-up focuses on the details of 
the form factor calculation and consequences; however, we will comment on contributions to b —> s 
from nonlocal matrix elements and other open issues at the end. 

We used a subset of the MILC Collaboration gauge field configurations which include effects 
of 2 + 1 flavors of 0(a 2 ), tadpole-improved (asqtad) sea quarks. We chose three ensembles, al¬ 
lowing for modest tests of discretization and quark mass effects. We correctly anticipated that 
obtaining a good signal over statistical noise would be a challenge. In order to address this issue 
we used eight light quark sources on each configuration, obtaining over 30000 estimates of each 
correlation function on each of the three ensembles. 

The calculations were done with nonrelativistic b quarks, formulated in the Bt s \ rest frame 
and accurate through 0(v 4 ). The matching of the effective field theory currents to the physical 
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ones, done to ,GC s Aqct>/ mb, A.Q CD /m|), is the source of the largest controlled systematic un¬ 

certainty. 

The momentum of the K* or 0 meson varied from 0 to no greater than 4k/ L (~ I GeV) in 
magnitude, primarily due to signal-to-noise degradation. Thus the lattice calculations are done in 
the low recoil regime, i.e. at large q 2 , the lepton-pair invariant mass-squared. 

From the imaginary-time correlation functions, we performed bootstrapped tits to extract nu¬ 
merical estimates of the seven linearly independent form factors. The basis we use is {V.Ao, A\, 
A 12 , 7j, 7?, T 23 } (see [1] for definitions). Once these are determined on each ensemble, for sev¬ 
eral values of final state momentum, it remains for us to parametrize the shape and give results 
corresponding to the physical limit. 


3. Updated fits 


This section describes a minor update of our published results [1]. We fit each form factor, 
generically denoted F(t), using the following parametrization 


F (t) = 1 - t /(mg +A m F ¥ 0 + c ()\^ x + 4\s^s) +ofz(f;f 0 )] (3-1) 


where Ar = (m 2 - >n 2 n phy J/{47if K ) 2 and Ax, = {m\ - »^ s , p h ys )/( 4 xfn) 2 . Thus the parameters 
Cq| and Cq| v quantify the dependence of the form factors on the light and strange quark masses, 
respectively. The splitting Am F between the B,- Sj mass and the relevant resonance is taken to be a 
fixed value, but variations of up to 20% had no effect on the resulting physical form factors. We 
determine Cq 1v from a combined lit to B -y K *, B s —y 0, and B s —y K* data [1] and include it as a 
Gaussian prior in the subsequent fits to determine a{(, a F , and c r ()] . Physical results are obtained 
by using our fit results for a/- and a F with At = Ar v = 0. We tried several fits which allowed for 
discretization effects, but no statistical signal was seen for the corresponding parameters. 

At two kinematic points the form factors are not all linearly independent. When q 2 = 0, 
equations of motion may be used to relate matrix elements, giving two constraints on the form 
factors. At the kinematic endpoint, q 2 = t- = {m B — my) 2 , two pairs of form factors become 
linearly dependent: 


^12(0) 


2 2 

m B ~ m V 

Sm B my 


-A„(0) 


7i(o) = r 2 (o) 




723(G) 


(m B +m v )(m 2 B -ml-t-) 

16m Bitty 1 

(m B + m v )(ml + 3 m\ -1-) 

Sm B m 2 2 ^ 


(3.2) 


In our published work fits were done to individual form factors separately, with the exception of 
a joint fit to 7) and 7), so that the constraint 7) (0) = 7F (0) could be implemented. This was the 
most important constraint for existing measurements such as 38(B —> K*y). In these proceedings 
we implement all constraints in (3.2) by performing simultaneous fits to the vector and axial vector 
form factors {V,Aq,Ai,Ai 2 } and to the tensor and pseudotensor form factors {T\ This 

will ensure that endpoint relations such as F B = i are precisely satisfied [13]. The constraints are 
implemented by adding a fake “data point” to the fit requiring that the left-hand and right-hand 
sides (in the physical mass limit) are equal up to some uncertainty. For q 2 = 0 this uncertainty is 
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Figure 1: Lattice QCD determinations of B —»• K* (left), B s —>• (j) (middle), and B s —>• K* (right) form factors 
in the physical limit. Error bands include statistical and systematic uncertainties. 

set to be 1%, an estimate for the size of Oiz 2 ) terms truncated from our fit ansatz. For q 2 = q 2 wx 
we take the uncertainty to be 10 4 since we have data in the high -q 2 range. In practice, results do 
not change if either of these is reduced. 

We have tried including cross-correlations between all of the form factor data in order to 
estimate the correlations between the full set of fit parameters. The data covariance matrix was not 
determined well enough to allow for a single fit, so several fits to sets of 3 or 4 form factors were 
done. We suspect that the data covariance matrices were still not entirely well-determined. For 
example, when we included cross-correlations but no new kinematic constraints compared to our 
published fits, we found deviations of approximately 1 standard deviation, despite the fact that the 
fit parameters were independent of each other. Therefore, for our final fit we include only the data 
correlations between Aq and Au and between 7) and 73. These are the pairs of form factors which 
must satisfy kinematic constraints (3.2) at q 2 = 0. Figure 1 shows the results for the form factors 
in the physical limit. Further details and results of these fits are appended to the end of the arXiv 
version of this write-up [arXiv: 1501.00367], 

We use these new form factor fits, including the constraints (3.2) and the estimated correlation 
matrix, to determine several B {) -> X'°/i /J , B ± —> K + j .1 /j. , and B s —> 0p /J observables. 
For most quantities, notably the differential branching fraction, we find only negligible changes 
to our published results fits [2], Fig. 2 shows our Standard Model differential branching fraction 
compared to experimental data, with the data consistently lower than the theoretical prediction (see 
also Table 1). Tables 2 and 3 give Standard Model predictions for angular observables. The most 
significant change compared to [2] is to F/ , which shifts and is more precisely determined. This 
is due to the inclusion of the A^itJ) /A\{t-) constraint. There is also a < la shift in the central 
value of S 3 in the case of B s -> 0/x + /i~. 

We have repeated our beyond-the-Standard-Model fit to large -# 2 experimental data for 7?° -y 
1 /J and B s —> allowing the Wilson coefficients Cg = Cf M + C^ P and C[ } to deviate 
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Figure 2: Differential branching fractions for B° —> K*°ji 1 fi (left), B ± —> fC‘ + fi 1 /i (middle), and B s —)• 
$11+fl (right). Experimental results shown are from [3, 6, 5], respectively. 


q 2 bin (GeV 2 ) 

B° -> K*°ii+ii- 

B ± k*+ ii+ii~ 

B s (j) fi+ll ~ 

[14.18,16.00] 

[16.00,19.00] 

[14.18,19.00] 

7.8(1.2) x 10~ 8 
5.73(79)x 10~ 8 
6.52(94) x 10~ 8 

8 .4(1.3) x 10~ 8 
6.19(85)x 10~ 8 
7.0(1.0) x 10" 8 

7.7(1.0) x 10~ 8 
5.14(67)x 10” 8 
6.11(80)x 10~ 8 


Table 1: Standard model predictions for branching fractions d3S/dq 2 {GeV 2 ) in bins of q 2 . 


q 2 bin (GeV 2 ) 

F l 

Afb 

s 3 

P> 

M 

P> 

r 5 

[14.18,16.00] 
[16.00,19.00] 
[14.18,19.00] 

0.360(42) 

0.336(25) 

0.347(32) 

0.410(35) 

0.347(21) 

0.375(26) 

-0.160(29) 

-0.230(17) 

-0.198(22) 

0.612(17) 

0.650(08) 

0.633(12) 

-0.702(59) 

-0.541(35) 

-0.614(45) 


Table 2: Standard model predictions for B° —> K' {) u ' Li angular observables in bins of q 2 . 


q 2 bin (GeV 2 ) 

Fl 

s 3 

P’a 

[14.18,16.00] 

0.382(20) 

-0.172(13) 

0.624(8) 

[16.00,19.00] 

0.347(12) 

-0.242(08) 

0.659(4) 

[14.18,19.00] 

0.364(15) 

-0.209(10) 

0.642(5) 


Table 3: Standard model predictions for B s —> (till 1 ji angular observables in bins of q 2 . 

from their Standard Model values (Cg SM ~ 0). This yields C^ p = —1.1 ±0.5 and C 9 = 1.2 ±0.9, 
comparable to what we found before [2]. Figure 3 shows a contour plot of the likelihood function 
of this fit. 

4. Open issues 

One uncontrolled approximation in our calculation is that we fit correlation functions using a 
single interpolating operator for the vector meson final state, assuming that it corresponds to the 
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Figure 3: The likelihood function of a 2-parameter fit to B° —>• K' ( \u jl and B s —t (!>Li 1 u experimental 
data with q 2 > 14.18 GeV 2 . The Standard Model (SM) value lies just outside the 2cr contour. 

K* or 0. A fully controlled calculation would include scattering states in the analysis. In order 
to do so, a much more expensive and complicated set of calculations must be undertaken. The 
path forward was set out in a paper which appeared during this conference [14]. Work to study the 
spectrum of Kn and Krj states has recently begun [15, 16, 17], marking the first major step toward 
a full lattice calculation of B —> K*(—> Kn) matrix elements. In the interim, the use of our form 
factor results comes with the assumption that threshold effects will be small. One might expect 
them to be smaller for B s —> 0 than for B —y K* since the 0 is relatively narrow. One might also 
note that heavy meson chiral perturbation theory predicts percent-level threshold effects in B —> D* 
form factors [18, 19] and hope that such is the case for lighter-mass mesons. Nevertheless at least 
there is a plan to systematically include this effect in the future. 

Perhaps a larger open question regards the extent to which form factors are sufficient to deter¬ 
mine the hadronic contributions to observables. In addition to matrix elements of 2-quark operators 
(2.2), there are matrix elements of non-local operators Oi(x)j^(y), notably where (), is a 4-quark 
b —> s operator which creates a cc pair, annihilated at another point by the vector current / u . It was 
expected that these matrix elements would make theoretical predictions unreliable for \Jq 2 ~ mjy 
or my. Within the context of an operator product expansion, leading-order contributions from these 
non-local operators can be (and, in our work, have been) included [20, 21]. However, contributions 
from the i/r(4160) to the B —> K 1 /d ‘ /d~ are larger than anticipated [22] and the same appears to 
be true for decays to vector mesons. Understanding nonfactorizable contributions is probably the 
most important open issue which theoretical predictions must confront [23] . 

5. Conclusions 

This write-up summarizes our recent calculation of B —y K*, B^ -y 0, and B s —» K* form fac¬ 
tors. Here we give a slightly improved set of fits to the form factor shapes, including all four 
kinematic constraints and estimates of the more significant correlations between form factor pa¬ 
rameters. Given that this update is only a minor improvement upon our published work and has 
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a small effect on observables, we would be grateful if those authors using the form factor results 
presented here would cite Ref. [1] in addition to this work. 
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Form factor F 

B -y K* 

B s -A (f) 

B s -A K* 

Ao 

87 

0 

-87 

V, F 

135 

45 

-42 

Al, A 12 , T 2 , 723 

550 

440 

350 


Table 4: Mass differences Am 1 (in MeV), between the initial state and pertinent resonance (or effective 
pole) contributing to form factor F. 


A. Supplemental material 


In this Appendix we give detailed tables in order for readers to use our new fits to reconstruct 
the form factors (in the physical limit). We also compare the new fits to the published ones [1]. 

Table 4 gives the mass differences between the B or B s and the resonance (or effective pole) 
relevant for the various form factors. These enter the prefactor of the form factor shape (3.1). The 
remaining tables give results for the constant and linear coefficients (ao and a\) of z in the form 
factor shape (3.1): Tables 5 and 6 for B -y K* decays and Tables 7 and 8 for B s -a <j) decays. In 
Tables 9 and 10 we give results for B s —* K* decays, such as the b -y u decay B s -A- K*lv and the 
b —> d decay B s — > K*il. These fit results come from separate fits for each Table, with correlations 
included only between Ao and A 12 and between 7) and 73. Other data correlations were not included 
for reasons discussed in Sec. 3. 1 

Those wishing to use our form factor calculations for predictions of observables need not 
include all the correlation data. Large q 2 observables predominantly depend on V and A 1 , and 
their fit parameters are not very correlated. On the other hand, the parameters for individual form 
factors, £(o and a 1 , arc highly correlated, so the uncertainty of the LQCD-determined form factor 
F(t ) [Eq. (3.1)] is given (in the physical limit) by 


SF(t) = 


1 -t/(m B(s) +Am F ) 2 


(Safi) 2 + (8a\)~z{f,to) 


+ 


(5flo)(5flf)z(f;t 0 ) 


(A. 1) 


Since physical results are obtained by setting Ax = Ax s = 0 in (3.1), one does not need our results 
for coi or coi. s except to look at the quark-mass dependence of the form factors. Data files containing 
the fit results as tabulated in Tables 5-10 are included as ancillary files with the arXiv submission. 

Figures 4, 5, and 6 show the fit results presented here compared to our published results [1]. 
The most visible difference in the Figures is the reduction of uncertainty in Ao at low values of q 2 
(positive z) due to the constraint that it equals A 12 at q 2 = 0 up to a known multiplicative factor. 
The inclusion of the q 2 = A constraints ensures endpoint relations are accurately satisfied [13]. 
Together these improvements to the fits yield slightly more precise and accurate predictions for 
several observables. Table 11 gives values for each of the form factors at a few reference kinematic 
points. 


1 ArXiv version 1 of this write-up gave an estimate of the full correlation matrix using fits to subsets of form factors. 

However, the resulting matrix was not positive semi-definite. One could devise a method for removing the negative 
eigenvalues, but since the differences between the fits given in versions 1 and 2 of these proceedings are not significant 
compared to other uncertainties, we believe the neglected correlations are presently insignificant. 
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p 

value 

C p< 

C pA 

C A 0 
PA n 

C A 0 
P>«1 

C A, 

w 

C A1 C Ai 9 

< 

0.4975(667) 







a\ 

-2.015(916) 

0.8590 






A 0 

«0 

0.5023(370) 

0.0000 

0.0000 






-1.608(447) 

0.0000 

0.0000 

0.6673 




Ai 

«0 

0.2848(233) 

0.0000 

0.0000 

0.0002 

-0.0401 



af 1 

0.191(280) 

0.0000 

0.0000 

0.0000 

0.0104 

0.9483 


tf 4 ' 2 

“0 

0.2195(238) 

0.0000 

0.0000 

0.9043 

0.9091 

0.0073 

-0.0019 

An 

flj 

0.332(300) 

0.0000 

0.0000 

0.8878 

0.9272 

-0.0487 

0.0126 0.9756 


Table 5: Fit parameters and correlation matrix elements for B —t K* form factors V, Ao, Ai, and A\ 2 - 


P 

value 

C n 

pa 0 

C n 

PA 1 1 

C t 2 
PA n 

C r, 
PA 1 

C 793 

pK 5 

T, 

a 0 

0.4197(241) 






T, 

Cl\ 

-1.363(259) 

0.5005 





t 2 

a 0 

0.27997(1948) 

0.8503 

0.8015 




t 2 

a l 

0.117(236) 

0.8205 

0.8364 

0.9324 



T23 

a 0 

0.5235(451) 

0.0064 

-0.0022 

0.0117 

-0.0061 


T 22 , 

a x 15 

-0.271(579) 

-0.0341 

0.0118 

-0.0628 

0.0325 

0.9520 


Table 6: Fit parameters and correlation matrix elements for B —y K* form factors 7 ], If and 7 23 . 


P 

value 

C pA 

C pa\ 

C A 0 
PA 0 

C A 0 

PA , u 

C Aj 
PAf 

C Ai C Ai 9 

PV\ PA 0 

a V o 

0.4525(303) 







a\ 

-2.399(496) 

0.8044 






A 0 

V 

0.52935(1633) 

0.0000 

0.0000 





af° 

-1.592(243) 

0.0000 

0.0000 

0.7791 




Ai 

V 

0.28283(1107) 

0.0000 

0.0000 

-0.0097 

-0.0774 



1 

0.1238(1501) 

0.0000 

0.0000 

-0.0006 

-0.0047 

0.9353 


fl 4 * 2 

“0 

0.20661(978) 

0.0000 

0.0000 

0.9322 

0.9364 

0.0022 

0.0001 

flf 12 

0.4763(1473) 

0.0000 

0.0000 

0.9111 

0.9598 

-0.0904 

-0.0055 0.9690 


Table 7: Fit parameters and correlation matrix elements for B s —> @ form factors V, Ao, A1, and A12. 


P 

value 

C 7, 

PA n 

C 7, 

P.O| 

C t 2 
PAn 

C 7 , 

PA\ 

C r 23 

T, 

a 0 

0.40160(1057) 






T, 

a l 

-1.1340(1325) 

0.3041 





t 2 

a 0 

0.28297(752) 

0.7021 

0.6263 




fl[ 2 

0.2487(975) 

0.6767 

0.6467 

0.9273 



?23 

°o 

0.52193(1566) 

0.0202 

0.0104 

0.0351 

0.0080 


?23 

0.384(236) 

-0.0689 

-0.0357 

-0.1199 

-0.0273 

0.9455 


Table 8: Fit parameters and correlation matrix elements for B s —> @ form factors 71, 72 , and T 33 . 
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p 

value 

C p< 

C pA 

c A 0 
P<“ 0 

C A 0 
P,a, 

C A! 
P’° 0 

C Ai C A19 

P& 1 P0 0 

a V o 

0.3367(466) 







a\ 

-2.879(661) 

0.9037 






Aq 

a 0 

0.50344(1856) 

0.0000 

0.0000 





fl f° 

-1.855(235) 

0.0000 

0.0000 

0.5757 




Ai 

a o 

0.23390(1151) 

0.0000 

0.0000 

-0.0173 

-0.0681 



A ! 
a l 

0.0907(1316) 

0.0000 

0.0000 

0.0140 

0.0554 

0.9119 


A P 

«0 

0.20303(1233) 

0.0000 

0.0000 

0.8777 

0.8804 

-0.0175 

0.0142 

af 12 

0.4442(1509) 

0.0000 

0.0000 

0.8480 

0.9095 

-0.0785 

0.0638 0.9637 


Table 9: Fit parameters and correlation matrix elements for B s K form factors V, Ao, Ai, and A 12 . 


P 

value 

C T t 

P’ a 0 

C r. 

C r 2 
P' a n 

C r. 

C t 23 

P> a 0 

Ti 

a o 

0.34878(1423) 






a\ l 

-0.9807(1783) 

0.1851 





Tt 

a 0 

0.24243(1063) 

0.6786 

0.6929 




Ti 

a l~ 

0.2051(1220) 

0.6440 

0.7236 

0.9296 



T23 

a 0 

0.4701(237) 

-0.0026 

0.0018 

-0.0052 

0.0044 


T23 

0.138(291) 

-0.0356 

0.0244 

-0.0703 

0.0595 

0.9468 


Table 10: Fit parameters and correlation matrix elements for B s — > K form factors T \, T 2 , and 'A 3 . 





B 

^K* 




9 2 

V 

Ao 

A, 

An 

Ti 

t 2 

Ti3 

7 max 

1.92(15) 

1.90(13) 

0.620(35) 

0.444(25) 

1.54(10) 

0.623(35) 

1.255(71) 

16 

1.28(11) 

1.280(90) 

0.523(36) 

0.389(30) 

1.049(70) 

0.520(33) 

1.010(67) 

12 

0.84(12) 

0.861(77) 

0.440(42) 

0.339(40) 

0.711(54) 

0.433(37) 

0.809(81) 

0 

0.31(15) 

0.351(74) 

0.303(51) 

0.251(53) 

0.291(44) 

0.291(44) 

0.50(10) 

<r 

V 

Ao 

B 

A, 

s ^ 0 

A12 

Ti 

t 2 

T 23 

7 max 

1.74(10) 

1.856(98) 

0.624(32) 

0.396(21) 

1.351(75) 

0.605(31) 

1.128(58) 

16 

1.195(73) 

1.325(71) 

0.529(29) 

0.359(20) 

0.979(54) 

0.521(27) 

0.966(51) 

12 

0.767(64) 

0.907(53) 

0.439(28) 

0.321(22) 

0.680(38) 

0.439(25) 

0.810(47) 

0 

0.244(71) 

0.391(40) 

0.294(28) 

0.248(26) 

0.303(23) 

0.305(22) 

0.555(45) 

7 2 

V 

Ao 

B s 

A, 

-a-K* 

An 

Ti 

t 2 

Ti3 

7max 

1.99(13) 

2.35(13) 

0.582(32) 

0.424(23) 

1.470(98) 

0.578(31) 

1.180(63) 

16 

1.035(78) 

1.355(77) 

0.451(26) 

0.363(22) 

0.890(56) 

0.459(26) 

0.910(52) 

12 

0.584(86) 

0.884(55) 

0.370(26) 

0.321(25) 

0.605(39) 

0.383(25) 

0.743(53) 

0 

0.07(10) 

0.335(39) 

0.242(26) 

0.243(29) 

0.259(27) 

0.261(25) 

0.483(55) 


Table 11: Form factor values (and total estimated uncertainties) at several reference values of q 2 (in GeV 2 ) 
using updated fits to our lattice results. 
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Rare B decays using lattice QCD form factors 


M. Wingate 




zlt, 12 GeV 2 ) z(t, 12 GeV 2 ) 



zlt, 12 GeV 2 ) 



zlt, 12 GeV 2 ) 




z(t, 12 GeV 2 ) zlt, 12 GeV 2 ) 


Figure 4: Comparison of B —> K* form factor fits presented here (blue dashed and dotted lines) compared to 
published fits [1] (teal bands). Black symbols are LQCD results at unphysical quark masses, whereas curves 
and bands are extrapolated to the physical limit. 
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Rare B decays using lattice QCD form factors 


M. Wingate 
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z(t, 12 GeV 2 ) 
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z{t, 12 GeV 2 ) 



0.00 
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Figure 5: Comparison of B s —> 0 form factor fits presented here (blue dashed and dotted lines) to published 
fits [1] (pink bands). Black symbols are LQCD results at unphysical quark masses, whereas curves and 
bands are extrapolated to the physical limit. 
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Rare B decays using lattice QCD form factors 


M. Wingate 



zft, 12 GeV 2 ) 



2 ft, 12 GeV 2 ) 



2 ft, 12 GeV 2 ) 



2 (t, 12 GeV 2 ) 



2 (t, 12 GeV 2 ) 



2 ft, 12 GeV 2 ) 


Figure 6: Comparison of B s —> K* form factor fits presented here (blue dashed and dotted lines) compared 
to published fits [1] (orange bands). Black symbols are LQCD results at unphysical quark masses, whereas 
curves and bands are extrapolated to the physical limit. 
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